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Abstract. Based on the Gaussian wave packet solution for the harmonic oscillator and the 
corresponding creation and annihilation operators, a generalization is presented that also 
applies for wave packets with time-dependent width as they occur for systems with different 
initial conditions, time-dependent frequency or in contact with a dissipative environment. 
In all these cases the corresponding coherent states, position and momentum uncertainties 
and quantum mechanical energy contributions can be obtained in the same form if the cre- 
ation and annihilation operators are expressed in terms of a complex variable that fulfills a 
nonlinear Riccati equation which determines the time-evolution of the wave packet width. 
The solutions of this Riccati equation depend on the physical system under consideration 
and on the (complex) initial conditions and have close formal similarities with general su- 
perpotentials leading to isospectral potentials in supersymmetric quantum mechanics. The 
definition of the generalized creation and annihilation operator is also in agreement with a 
factorization of the operator corresponding to the Ermakov invariant that exists in all cases 
considered. 



1 Introduction 

Shortly after the first communication on wave mechanics pQ Schrodinger himself already 
considered the possibility of constructing stable wave packet (WP)-type solutions of his 
equation that behave closest to the classical particles in the case of the harmonic oscillator 
(HO) [2j. The functions that fulfill this requirement, the Gaussian WPs, are determined 
by two parameters, the maximum and the width, where in his case the maximum followed 
the classical trajectory and the width was just a constant. This led him to think that the 
WP could just describe a distribution of matter instead of a mass point, but this hope 
was destroyed by the rapidly spreading WP for the free motion. 
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Generalizations of Schrodinger's approach were achieved in the description of coherent 
light beams emitted by lasers [3H5], considering what is now known as coherent state 
(CS). There are at least three different formal definitions of these states in the literature 
[6]: (1) minimum uncertainty CS, meaning Gaussian WPs that minimize the Heisenberg 
uncertainty relation; (2) annihilation operator CS and (3) displacement operator CS. The 
definition of the annihilation operator CSs involves the factorization of the Hamiltonian 
(of the HO) in terms of creation and annihilation operators, a procedure that was also 
applied by Schrodinger [7], but already mentioned in principle by Dirac [8]. All this works 
perfectly well for at most quadratic and conserved (time-independent) Hamiltonians with 
Gaussian WP solutions. However, this type of WP solutions also exists if the Hamiltonian 
is no longer a constant of motion, but explicitly time-dependent (TD), e.g. like in the 
case of a HO with TD frequency u(t). In this dynamical invariant (with the 

dimensions of an action instead of energy) still exists, as has been shown in a formal 
mathematical context by Ermakov [9] already 45 years before quantum mechanics was 
developed by Schrodinger and Heisenberg. This invariant can be obtained by eliminating 
the TD frequency u(t) from the Newton equation of motion for the WP maximum with the 
help of an auxiliary equation that, as will be shown later, determines the time-evolution 
of the position uncertainty (or WP width). The important role of this Ermakov invariant, 
particularly in a quantum mechanical context, has been pointed out in the works of Lewis 
and Riesenfeld [10] and has been exploited by many other authors. In this context, also 
generalized CSs for these systems have been discussed [TT] . 

In this work we will give a further generalization that not only covers Hamiltonians 
with TD (harmonic) potentials, but also TD Hamiltonians that describe dissipative sys- 
tems where the energy is no longer a constant of motion and other, particularly nonlinear 
(NL) extensions of the TDSE, that try to include the effect of a dissipative environment 
into a quantum mechanical context. 

It is assumed that also for these systems a description in terms of pure states is possible 
where these pure states fulfill modified SEs (with possibly non-Hermitian Hamiltonians) 
that still have exact Gaussian-type WP solutions. There are several of these approaches in 
the literature [T2TH6] and we will show how some of the most accepted ones are interrelated 
and the meaning of their physical interpretation. For some of these approaches, it can be 
shown [T7] that also in these cases an exact dynamical invariant of Ermakov-type exists. 
This invariant can be expressed in different variables corresponding to the respective 
Hamiltonian, but if the proper transformations between these approaches are taken into 
account, they all lead to the same invariant. 

This invariant, as in Schrodinger's original approach and in the extension to include 
TD frequencies for the HO, can also be factorized and, for the quantized version of the 
invariant, this leads to generalized creation and annihilation operators that are also ap- 
plicable for open dissipative quantum systems. 

A key point of our approach is that we do not express these generalized creation 
and annihilation operators in terms of the (real) auxiliary variable that fulfills the (NL) 
Ermakov equation, but in terms of a complex variable that fulfills a quadratic NL Riccati 
equation. The imaginary part of this variable satisfying the Riccati equation is uniquely 
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related to the position uncertainty or the real Ermakov variable. The particular form of the 
real part can change according to the system of interest (conservative or dissipative) . In 
our generalized formalism, the complex Riccati equation exists (at least) for HO systems 
with constant or TD frequency and with or without linear velocity dependent friction 
forces (where the damped free motion is covered for oj — > 0). So, the same form of the 
generalized creation/annihilation operators in terms of the Riccati variable is applicable 
for all these systems which is not the case if they are expressed in terms of the real 
Ermakov variable. 

A further advantage of the complex Riccati approach over the real Ermakov one is 
that in the Riccati case, due to the quadratic nonlinearity (particularly in the dissipative 
case), the possibility of pairs of solutions with different physical properties (connected 
with phenomena like Hopf bifurcations) immediately becomes obvious, which is not at 
all the case if the corresponding real NL Ermakov equation is considered. The solutions 
of the Riccati equation are quite sensitive to the (complex) initial conditions and can, 
for even slightly different ones, lead to a qualitatively very different behaviour of the 
quantum system. Due to the Riccati equation, formal similarities with the supersymmetric 
approach to quantum mechanics [T5H2U] can also be established and exploited. 

We introduce our generalized CSs in terms of the complex Riccati variable using the 
TDSE for the HO with the possibly TD frequency as starting point and show, that our 
approach is consistent with the definitions of displacement operator CSs and annihilation 
operator CSs. Since we consider general CSs, with position-momentum correlations, the 
width of the Gaussian WP solutions is no longer constant and these CSs do not fulfill the 
minimum uncertainty requirement, but the modified one associated with the Schrodinger- 
Robertson uncertainty relation [2T] . 

The approaches to include dissipation and irreversibility that will be discussed have one 
thing in common: they lead to Langevin-type evolution equations (with linear velocity 
dependent friction force) for the classical part of the system, determining the motion 
of the WP maximum. No stochastic force appears explicitly, but this aspect might be 
taken into account by non-Hermitian imaginary contributions to the Hamiltonian in the 
NLSE approaches. Two of these approaches start already on the quantum level and 
add some (NL) terms to incorporate the effect of the dissipative environment (leaving 
the standard definitions and interpretations of operators and wave functions unchanged). 
Two other ones modify the canonical Lagrange/Hamilton formalisms on the classical 
level in a way that allows one to include the above-mentioned friction force, but after 
quantization caution must be taken with respect to the physical meaning of operators and 
wave functions in this canonical context. We will show that all the considered approaches 
share the same Ermakov invariant and how this can be expressed and factorized in physical 
terms. 

In particular, the physically-relevant Riccati equation, its possible solutions and rela- 
tions to supersymmetric quantum mechanics will be discussed and the quantum uncer- 
tainties will be expressed in terms of the complex Riccati variable. 
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2 Gaussian wave packets, complex Riccati equations 
and generalized coherent states 

2.1 Gaussian wave packets and Ermakov invariant 

In the following, TDSEs with at most quadratic Hamiltonian (particularly, the HO with 
constant or TD frequency u, and the free motion, V — in the limit u — > 0) in one 
dimension shall be considered, 

d f ft 2 d 2 171 1 

In these cases, Gaussian wave packet solutions can be obtained that can be written in the 
form 

1 



*wpOM) = iV(t)exp <^ % 



y(t)x 2 + -(p}x + K(t) 



(2) 



with x = x— (x) = x—i] where the mean value of position is given by (x) = J*_^°° dx^/*x^/ = 
rj(t), (p) = mi] and the purely TD function K(t) and normalization factor N(t) (which will 
be specified later). The time-dependent coefficient of the quadratic term in the exponent 
is assumed to be complex, y(t) = Ur + iyi where the imaginary part is related to the 
position uncertainty (x 2 ) = (x 2 ) — (x) 2 via yj = ^^y- The maximum of the WP is 
located at x — (x)(t) = rj and, thus, follows the classical trajectory determined by the 
Newtonian equation of motion that can be obtained by inserting (2) into (1), 

ii + u 2 (t)r] = 0. (3) 

The time-evolution of the WP width can be determined in the same way and is governed 
by the complex nonlinear Riccati equation 

^VUW+^O. (4) 



m J \m 

So, solving the classical equation of motion (3) and equation (4) provides all the informa- 
tion that determines the time-evolution of the quantum mechanical WP (2). 

The complex variable of the Riccati equation will play the key role in our 

construction of the CSs. Since, for the HO with TD frequency and the dissipative systems 
the Hamiltonian is no longer a constant of motion, a different quantity shall be used for 
the factorization method; namely, the Ermakov invariant. This invariant can be obtained 
by eliminating the frequency u from Eqs.(3) and (4). Usually, for this purpose, Eq.(4) is 
written in the form of a (real, NL) Ermakov equation that can be obtained by introducing 
a new variable a(t) via 

2h 1 M 
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where a = \J^y{x 2 ) is directly proportional to the WP width. Inserting this definition 
(5) into the imaginary part of Eq.(4) yields 



2h a 
-Vr 



(6) 



m a 

and, finally, inserting (5) and (6) into the real part of (4) leads to the Ermakov equation 



a + u 2 (t)a = —5-. 

cr 



(7) 



Via the afore-mentioned elimination process, the invariant can be obtained as 

(fi(\ - nn)~ ----- I — ! = — n" i /• - — ;/ ! 

1 



(rja - ar]) 2 + 



= 2 a 



rj 



2h \ V ((ih . 

— Vr )v) + — Vi V 
m II \\m 



(8) 



const. 



The corresponding quantum mechanical operator can be obtained easily by replacing f\ 
by the momentum operator p op — | divided by m and r\ by x, taking into account that 
p op and x do not commute. 



2.2 Complex Riccati equation and supersymmetric quantum 
mechanics 

There are different ways of treating the Riccati equation (4); in the following, we will 
use a method that best shows the similarity to supersymmetric quantum mechanics and 
point out the properties that have qualitative consequences for the corresponding WP 
solutions. For this purpose we take advantage of the fact that the inhomogeneous Riccati 
equation can be transformed into a homogeneous Bernoulli equation if a particular solution 
(m^) °f ^ ne R-i cca ti equation is known. The general solution of Eq.(4) is then given by 



2h 



y = ^y+ ™ v (t) where ^v(t) fulfils 



2h.\ {AhA f2h \ f2h x 2 

—v + — y —v + —v 
m J \m J \m J \m 



0. 



(9) 



The coefficient A = 2 (^y) of the linear term depends on the particular solution; further 
details will be discussed later. The NL Bernoulli equation (9) can be linearized via 



m wyt) 



to yield 



w — Aw = 1, 



which has the solution 



w(t) 



w + 



r* 

dt'e-f 



dt"A 



j'dt'A 



(10) 
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Particularly for constant A, this reduces to 

w(t) = j{e At -l)+w e At , (12) 

depending on the (complex) initial value wq = (f^o) ■ Because Eqs.(4) and (9) are 
nonlinear equations, this dependence on the initial conditions can be very sensitive and 
lead to qualitatively quite different behaviour of the WP width. If, e.g., for the HO with 
constant frequency u = ojq the particular solution is chosen to be = i—yj = iuj 
(a solution with a minus sign is mathematically also possible but would results in an 

— 1/2 

unphysical WP solution), i.e., A = i2uo, a WP with constant width ao = u = 
(y^Io) 1 2 is obtained. Whereas, for any choice a ^ uj corresponding to — Vq = 
i _ ( i 



Uq) 7^ 0, an oscillating width is obtained (see also [22J). 



The similarity to supersymmetric quantum mechanics becomes obvious when we in- 
troduce the abbreviated form I{t) = J dt'e~J dt"A(t") ^ a }} ow j n g f or rewriting of the 
general solution of Eq.(4) as 

2h , . 2h^ d , r rr . 

—y{t) = —y + -\n[w + It, 13 
m m at 

defining a one-parameter family of solutions depending on the initial value of Wq = 
(|~Uo) as parameter. Comparison with supersymmetric quantum mechanics [T814201 
[23l[23] shows that solution (13) is formally identical to the most general superpotential 

W(x) = W(x) + —]n[X 1 +Xi(x)l (14) 
leading to a one-parameter family of complex isospectral potentials 

V^x) = W 2 - ^-W = V l - 2-^lii[A a + Z 1 (x)], (15) 

that have the same supersymmetric partner potential V^x) (see, e.g., [24TT29] ). In this 
case, the integral Ii(x) is defined as 

Ii(x) = / dx'^\{x) (16) 



where ^\{x) is the normalized ground state wave function of the SE with potential V\(x) = 
W 2 (x) — 4:W(x) and Ai is a (usually real) parameter. 

A major difference between this supersymmetric situation and the one in our case 
(apart from replacing the spatial variable by a temporal one) is the fact that the variables 
of the nonlinear Eqs.(4) and (9) are complex, whereas W(x),W(x) and Ii(x) are real 
quantities; also the parameter u> in our case is generally complex. This provides a larger 
variety but, nevertheless, certain methods and results can be transferred from one system 
to the other (which will be discussed elsewhere). 
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2.3 Generalized creation and annihilation operators and corre- 
sponding coherent states 

In the following it will be shown how the usual creation/annihilation operators for the HO 
with constant frequency Uq are related to the complex Riccati equation (4), how they can 
be generalized to take into account the additional dynamical aspects concerning the WP 
width and what consequences this has for the corresponding generalized coherent states. 

In the usual treatment of the HO in terms of creation and annihilation operators, the 
Hamiltonian 

1 m ( 1\ 

H = + = U + a + - , (17) 



with — |^ is divided by Tiwq to yield 



H= (a + a+- ) , (18) 



tvjjQ \ 2 

where H is dimensionless (but — = a^H would have the dimension of an action) and the 



creation and annihilation operators a + and a are defined as 

— - tUJ X ) 

. (19) 

+ • / m [Pop , ■ 

V 2nio V m 

The corresponding CS is a WP with constant width, corresponding to the particular 
solution of the Riccati equation with = i^yi = i^Vi = i^o- Replacing u by 4y or 
^Vi, respectively, the operators a and a + can be rewritten as 

— a ( Pop j( 2 ^ \ x 
2h I m \m J 

\\ (20) 

a+ =- i )j2h ao {m +i {m Vl]X 

As we have shown, already for u = uq = const., solutions with TD WP width exist, i.e., 
«o turns into a(t) and a ^ 0, hence —yR = -, must be taken into account. Obviously, 
the same also applies for the HO with TD frequency u = u(t). Therefore, in Eqs.(20 a,b) 
a must be replaced by a(t) and % (^yi) by (ffy), thus leading to 





2h \ m \ m 



3+( t ) = -i./^a (— ~ i—y" I * 



2h \m \m' 



(21) 
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It is easy to check that they satisfy the standard commutation relations if Eq.(5) is fulfilled. 

Since, at least for TD frequency u, the corresponding Hamiltonian is no longer a 
constant of motion, one might ask if a(t) and a + (t), as defined above, are constants of 
motion, i.e., if they fulfill 

|a + i|5,fl]. = (22) 

(with [ , ]_ = commutator). To answer this question, x, p and hence H shall be expressed 
in terms of a(t) and a + (t), leading to 



x 



V 



with 



rr ft 2 

H = —a 



(2h 

—VR 
\ m 



This finally leads to 



a \2fi 

— \ — (a + a 

2 \ m 



a 2h f (2h , 
m-\ — < —y 
I V m {\m 



'2h 

la + a T ) + % I — yi 

m 



d ~ 1 r~ rn 



a + — y a 
\ m 



(23) 



1 2 



-co a 



1 



a' 



a ± 0. 



(a + + a) . (24) 



(25) 



So, a(t) and a + (t) are no constants of motion but can be turned into such by simply 
introducing a phase factor according to 



r* mi J_ 



a(t) =a{t)e iI dt 
a+(t) =a + (t)e- i $ tdt 



t Ml 1 



(26) 



In order to elucidate the meaning of this phase factor we return to the Riccati equation 
(4) and apply another treatment, namely, the linearization using the ansatz 



V m 



A 
A' 



(27) 



u + iv 



ae 



%{p yielding the linear complex Newtonian 
= (28) 



with the complex quantity A 
equation 

X + u 2 (t)X 

for A(t), if (27) is inserted into Eq.(4). 

It should be mentioned that generalized creation and annihilation operators for the 
HO with TD frequency using a complex variable corresponding to A had been used by 
Malkin et al. [30]. Real and imaginary parts of A and their time-derivatives are also the 
elements of the 2x2 matrices of the real symplectic group Sp(2, R) that provides the 
representation of canonical transformations in TD quantum mechanics [31 j . 
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Using the polar form of A and inserting (27) into the Riccati equation (4), one obtains 
from the imaginary part a relation between the phase and amplitude of A as 



a- 



(29) 



So, the phase factor f di! \ is (up to a constant) just the angle (p(t) in the complex plane 
of the quantity X(t) that allows for the linearization of the Riccati equation (4). (For 
further details, see [221132].) 

Next, it shall be shown that one can use the definition of the generalized annihilation 
operator a(t) (Eq. 21a or 26a; the phase factor shall be omitted in the following since it 
can be absorbed in the purely TD function K(t) in the exponent of the WP/CS) or into 
N(t)) to create a CS \z) that is an eigenstate of a{t) with complex eigenvalue z, i.e., 



a(t)\z) = z\z). 



(30) 



First, the complex eigenvalue z shall be determined in terms of (x) z = r\ and (p) z = mi], 
where Eq.(30) is assumed to be valid. From 



—a{z + z) = \ —az R = r), 
2m V m 



h 

2m 



-am 



—y)z*+[ —y* z 
m J \m 



mfj, 



(31) 



follows that 




1 

7! 



2h 

-rot I —yj ] rj, 
a \m 



m 1 m 

-(^-r / d) = -^/-a 



2h 

—Vr } V 

m 



(32) 



or 




2h 



— )+ i (rja — rjdi) 
a 



Im 

v/2V¥ a 



2h \ ( (2h 
—yi h + l h _ —vr ) v 

ml \ \m 



(33) 



This shows the connection between the eigenvalues z (or z*) and the Ermakov invariant 



as 



/ 2 i 2 \ ^ * 

— [Zr + z R ) = — zz 

m v ' m 



m 



(34) 



Immediately from this it follows that the operator corresponding to the Ermakov invariant, 
when p = mx is replaced by the operator p op = |^ (and taking into account that 
[p, x\- = i), can be written in terms of our generalized creation and annihilation operators 
as Il, op = — [a + (t)a(t) + i] , which is in agreement with the approach of Hartley and Ray 
pTj to construct these operators for the HO with TD frequency from I L op via factorization. 
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In the position-space representation, the CS that is eigenstate of a(t), 

(x\a(t)\z) = z(x\z) or z^ z (x) = i\ ^a \ — - ( — y ) x \^ z (x), (35) 

V la [mi ox \m J J 

can be given as 

* g (x,t) = M(t)exp |^ (x - (x)) 2 + ~{p)x + 7^(p)(x) + \{z 2 + W 2 )} , (36) 

which is, for M(t) = A^(t)e"K^+l z l 2 ) and N(t) = (f|) 1/4 Q) 1/2 , identical to the normal- 
ized Gaussian WP (2), 

/ mx 1 / 4 / 1\ 1 ^ 2 f [ 1 1 ~n 

*z(x,t)={—) (-J ex V U yx 2 + -(p)x+—(p)(x) \ = V W p(x,t), (37) 

where, for the choice of a complex normalization factor N(t) = (fl) 1 ^ 4 (t) 1 2 j the TD 
function in the exponent has the form K(t) = ^(p)(x). Using the polar form of X(t), 

i.e., A = ae^ with if = fdt'^, N(t) can be written as N(t) = (^2) 1/4 e"^ /2 = 

{^Wm) 1/4e ~ qftdt '^- For « 2 = «o = "o\ N(t) turns into N(t) = (^f ) 1/4 e— °*/ 2 

and contributes the ground state energy of the HO that corresponds to E — ^(p 2 ) + 
yOJq(x 2 ) = |co> in the case of the WP with constant width. Here we obtain a general- 
ization for a = a(t). The phase factors e ±% $ dt ^ occurring in the creation/annihilation 
operators as defined in Eq. (26) can be absorbed into N(t) in a similar way. 

As mentioned in the introduction, the CS can also be defined as displaced vacuum 
state, i.e., 

\z) = exp {za + (t) - z*a(t)} |0) = D(z)\0). (38) 
In the position-space representation, the vacuum state 0o can be obtained via 

< x ^ = i M a {^-&) x } Mx) = (39) 

as 

(x\0) = <f> {x, t) = N(t)e iy{t)x2 . (40) 

Note that due to y(t) the exponent is now complex, in particular via iynx 2 , the term 
f £ 2 already occurs here naturally, whereas it must be introduced in the approach of 
Hartley and Ray [11] via a unitary transformation. 

Using the Baker-Campbell-Hausdorff formula, the CS can be written in the usual 
form 

nl 



\z) = D(z)\0) = e-sW exp(^a + (t))|0) = ^ 1 " ' |0), (41) 



n=0 
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now with a + {t) as denned in (21b) (or (26b)). Evaluating e za+ with the help of z and \z\ 2 , 
as given in (33) and (34), leads again in the position-space representation to ty z (x,t) as 
given in (37). 

Finally, it will be shown that our CS also fulfills the Schrodinger-Robertson uncer- 
tainty relation and how the uncertainties can be expressed in terms of yu and yi (since 
this form will also remain valid when dissipative effects are included, as will be explained 
in the next section). 

In terms of z, z* and the mean value of (x 2 ) can be written as 



(3*),= ± a 2 (z* 2 + z 2 + 2\z\ 2 + l) 
2m 

leading (together with Eq.(31a)) to the mean square deviation of position 

h 



(x 2 ) z = (x 2 ) z - (x) 2 z = ^a 2 . 



In the same way, from 



(2h 



I — y J z * 2 + I — y 

\ m I \m 



2h 



z 2 + 



2h 



m 



{2\z\ 2 + l) 



and using Eq.(31b), one obtains 



(42) 



(43) 



(44) 



(P 2 ), = (P 2 )* - (p)l = ^a 2 



2h 


2 


hm 






—y 




~2~ 




(-») 


m 






\ m J 



and from 

{[x,p] + ) z = (xp + px) z = ha 2 



(2h \ *n[2H \ (2h . 

— V }z +[—y)z+[ —y R 2z + 1) 
\m J \m J \ m 



with (31a) and (31b) the correlation uncertainty 
1 



From 



: \ 2 l X >P\ + 



h 2 



-a 



(2h \ 2 (2h 



— Vr + I —Vi 
\m I \m 



(45) 



(46) 



[X,P\ + ^ =(^[x,p] + ^ - (X)z(p)z = \o? (j^Vr) ■ ( 47 ) 



h 2 4 {2h 
~ T a Vr 

4 \m 



h 2 



(48) 



follows that our CS also fulfills the Schrodinger-Robertson minimum uncertainty condi- 
tion. 
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3 Complex Riccati equations, Ermakov invariants and 
coherent states for dissipative systems 

There exist several approaches for describing dissipative quantum systems using modified 
effective one-body SEs where the effect of the environment on the observable system is 
taken into account by dissipative friction terms without considering the individual degrees 
of freedom of the environment. This usually leads to NLSEs, SEs with explicitly TD 
Hamiltonians or non-Hermitian Hamiltonians. Several of these approaches are discussed 
in the literature [T2T[T6] . In this context it was also investigated [17] to ascertain if 
Ermakov invariants also exist for these approaches. There are two NLSEs [TSUTB] and two 
explicitly TD Hamiltonians [121133] for which this is the case and that will be discussed 
subsequently. It will also be shown that they are not independent of each other. The 
NLSEs are related in a way that is connected with a (unitary) phase transformation (for 
further details, see [31]; the NLSEs and the TD Hamiltonians, however, are connected via a 
non- unitary transformation [35J. Taking these facts into account, the Ermakov invariants 
of all four approaches are equivalent and can be transformed one into the other. 

In the following, these approaches will be presented whilst giving the corresponding 
WP solutions, equation of motion for the classical trajectory and the complex Riccati 
equation describing the respective time-evolution of the WP maximum and width, as well 
as the resulting Ermakov invariants. The corresponding creation/annihilation operators 
for the (physically relevant) NLSE will be defined and the CSs for the dissipative sys- 
tems obtained on this basis. Finally comparisons will be made with the case discussed 
previously. 

3.1 Explicitly time-dependent Hamiltonian of Caldirola and Kanai 

This approach [T2] starts on the classical level with an explicitly TD Lagrangian, 



L 



CK 



2 

With the corresponding canonical momentum 



-x 2 -V{x) e 7i . (49) 



d - 

P = -qtLck = mxe 11 = pe 7 *, (50) 

the Hamiltonian can be formulated as 

H C k = ^-e-^p 2 + e^V (x) , (51) 
2m 

which yields the proper equation of motion including a linear momentum (or velocity) 
dependent friction force with friction coefficient 7, 

d d 
p + 7P + -7T- V = rax + m-fx + — — V = 0. (52) 
ox ox 
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But, Hck neither represents the energy of the system, nor is it a constant of motion. 

It is important to realize that not only the canonical momentum p is different from 
the physical (kinetic) momentum p = met but, in particular, the transition from the 
physical variables x and p = mx to the canonical variables x = x, p = pe ft represents a 
non-canonical transformation. 

The transition to quantum mechanics is achieved by replacing, as usual, the canonical 
momentum with a differential operator according to p — > p op = leading to the 

corresponding Hamiltonian operator and, thus, to the modified SE 

d - * 

ih—^CK {X, t) = H CK ,op^CK (X, t) 



= <e-* 



(53) 

K 2 d 2 N 



2m dx 2 



+ e^V^ CK (x,t). 



For the systems considered in this work, this equation also possesses exact Gaussian WP 
solutions like (2). The equation of motion for the WP maximum is just one for the 
classical trajectory (including the friction force), in our notation 

i) + 7?) + to 2 i] = 0. (54) 

The modified Riccati equation for the complex variable ( 2 j}y) CK reads 



/2ft A ^ f f2h 



2 



— y) + e~ 7 * — y) + w 2 e 7 * = 0. (55) 

\ m J CK \ m J CK 

The imaginary part of this variable is again connected with the WP width via i^Vi) CK = 
2m(d?) CK ■> w here the subscript CK denotes the mean value being calculated with ^ck, the 
solution of Eq.(53). Introducing, like in the conservative case, a new (real) variable acx 
via {^yi) CK = i again allows for the transformation of this Riccati equation into a 
(real) Ermakov-type equation, 

e" 27 * 

acK + iacK + u 2 a C K = (56) 

a CK 

This equation, together with Eq.(56) for n, forms the required system of equations that 
possesses an exact Ermakov-type invariant, here given in the form 

2 / \ 2-i 



t C K-\ 



)*(: 



const. (57) 

acx J 

A major point of criticism concerning this approach is an apparent violation of the un- 
certainty principle. Defining the uncertainty product via 

Uck = (&)ck<?)ck = (x 2 )cK(p 2 )cKe-^ (58) 

by expressing the physical momentum p in terms of the canonical momentum p as p = 
e~ yt p, yields an exponential decay of this product. So, Uck can become smaller than ft 2 /4 
and even vanish for t — > oo. (A solution of this problem will be given later.) 
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3.2 Description in exponentially expanding coordinates 

This problem does not occur when applying a different non- canonical transformation [27] 
using an exponentially-expanding coordinate Q = e 7 */ 2 x and the corresponding canonical 
momentum P = mQ = me 7< / 2 (x + %x). 

The corresponding Hamiltonian function 

H exp = —P 2 + -n 2 Q 2 = const. (59) 

with Q 2 = (u 2 — is not only a constant of motion but, expressed in terms of the 
physical variables x and p as 



H exp 



1 2 7 m 1 2 

p H — px H uj x 

2m F 2 F 2 



= ±-pl + ^u 2 x 2 = E (60) 



shows that it also represents the initial energy E of the system (with p and xo being the 
initial momentum and position). 

From the Hamiltonian equations of motion one obtains the equivalent Newtonian 
equation 

Q + [uj 2 - Yj Q = (61) 

which, expressed in terms of the physical variable x, again provides the desired equation 
of motion including the friction force. 

Canonical quantization according to Q — » Q op — Q,P — > P op = f -§q provides the cor- 
responding SE with Gaussian WP solutions, only now the physical variable x is replaced 
by the canonical variable Q (and u 2 by Q 2 ). Apart from this, it looks exactly like the SE 
in the conservative case. 

Obviously, the corresponding Riccati and Ermakov equations can immediately be writ- 
ten as 

«tf) + (^V +tf .o (62) 

m J exp \ m J exp 

and 

a exp + tt 2 (t)a exp = -3—. (63) 

^exp 

Expressed in terms of a exp and Q, the corresponding Ermakov invariant can be written 
in the form 



i = l -a 2 

J-exp 2 ex P 



Q-—Q) 1 -[A-Q 



a 2 

^■exp 



const. (64) 



which is not only independent of u, but also does not explicitly contain the friction 
coefficient 7 (and therefore also exists for 7 = j(t)). 
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This approach is related to the afore-mentioned one by Caldirola and Kanai via a 
unitary transformation (on the classical level), Q = e 7 */ 2 x, P = pe" 7 */ 2 + m^xe 7 '/ 2 , and 
the Hamiltonians and action functions are connected via 

H exp = H CK + ^F 2 (65) 

with 

F 2 (x, P, t) = xPe 7 * /2 - m|x 2 e 7 ' (66) 

and 

S exp = S CK + m JxV* = (s + m|x 2 ) e 7 *. (67) 

Note that with Schrodinger's definition of the wave function \I/ via the action function, 
S c = |ln\l/, this canonical transformation corresponds to a (unitary) phase factor ^|x 2 
in the wave function. 

The quantized version of this approach causes no problem with the uncertainty prin- 
ciple [33] providing one consistently translates operators and wave function from the 
canonical level to the physical one. How this can be achieved becomes obvious after we 
present the approaches that already start on the quantum mechanical level and modify 
the Hamiltonian operator to include the effect of the dissipative environment, leading to 
NLSEs. 



3.3 NLSE with complex logarithmic nonlinearity 

The first attempt to introduce dissipation into the SE via a (real) logarithmic non- 
linearity by Kostin [T3] was based on obtaining the correct Newtonian equation in- 
cluding the friction force for the mean values. It however suffered several shortcom- 
ings [15j[T6], particularly, the reversible continuity equation for the probability density 
g(x,t) = ^* (x,t) \I/ (x,t) in the case of an irreversible dissipative system lacked reason- 
able explanation. Our approach [16] was therefore motivated by the attempt to break the 
time-reversal symmetry on all levels to include irreversibility, a phenomenon usually asso- 
ciated also with dissipative systems (but not necessarily occurring only in those systems). 
This can be achieved by introducing a diffusion term into the continuity equation thus 
turning it into a Fokker-Planck-type equation (in particular a Smoluchowski equation). 
Following a method by Madelung and Mrowka36, the (real) Smoluchowski equation can 
be separated into two complex equations, namely a (modified) SE for the wave function 
(x,t) and its complex conjugate, (x,t), if the separation condition 



Ko 

Dd^ = 1 ang-(\ng)) (68) 
Q 



with diffusion coefficient D is fulfilled. This leads to the NLSE 

ih—^ NL {x,t) = I + V{x)+ 1 - (ln^v L - (ln^vL» \ ^ NL (x,t) (G9) 
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with a complex logarithmic nonlinearity (for details, see e.g., [IE]). The real part of this 
nonlinearity depends on the phase of the WP, is identical to Kostin's term and influences 
the motion of the WP maximum. The imaginary part that corresponds to the diffusion 
term in the Smoluchowski equation has no influence on the WP maximum, but a strong 
one on the dynamics of the WP width. 

The equation of motion for the WP maximum is again identical to Eq.(54). The 
equation for the time-dependence of the width can now be obtained from the modified 
complex Riccati equation 



2h 



m 



NL 



'2h 
7 I —y 

m 



NL 



2h 



m 



+ u 2 (t) = 0. 



(70) 



NL 



Again, the imaginary part has the same relation to the WP width (or position uncer- 
tainty) as in the conservative case, (^Ui) NL = wpw • Introducing the variable a^L 



via Civr) 
form 



NL 



the real part of the complex Riccati-variable now takes the modified 



2h 

— Ur,nl 
m 



&NL 
OlNL 



7 

2' 



(71) 



With the help of these relations, the Riccati equation now turns into the Ermakov-type 
equation 

<V\l ■ \ ^" r I n vi = -5—, (72) 



a 



NL 



i.e., exactly the same equation as for a exp in the expanding system. Together with the 
dissipative equation for the WP maximum, Eq.(54), this leads to the exact Ermakov 
invariant 



l NL 



&NL 7 



&NL 



2 ]V 



a 



NL 



const. 



(73) 



3.4 General nonlinear "friction potential" 

A different attempt at finding a dissipative modification of the SE was used by Siissmann, 
Albrechtl4 and Hassel5 in trying to add possible combinations of position and momentum 
operators and their mean values to the linear SE in order to obtain a modification that 
provides the correct equation of motion (54) for the mean value of position, including 
the linear friction force. For this purpose they added some kind of "friction potential" 
W(x,p; \l/) to the linear SE, where W introduced a nonlinearity due to the occurrence of 
the mean values. 

The general form of this term is 

W Gen = 7 <p>(x - (x» + C 1 - [(x - (x», (p - (p»] + (74) 
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(with [ , ] + = anticommutator) that, for any choice of C, leads to the desired equation of 
motion for (x) = i](t). Particular choices for C lead to the following friction terms: 

Siissmann (C =1): Wsnu = % [(x - (x)),p] + (75) 
Albrecht (C =0): W Alb = j{p){x - (x)) (76) 

Hasse (C = 1/2) : W Has = *[(x- (x)), (p + (p»] + (77) 

Whereas the choices for C — 1 and C = lead to a wrong reduced frequency for the 
damped HO, C = \ provides the correct one. 

It should be mentioned that Whos can De written as a combination of the other two 
approaches as 

W HM = \ {W Siiss + W Alb ) . (78) 

The NLSE of Hasse also has Gaussian WP solutions for the systems under considera- 
tion and the equations of motion for its maximum and width (and therefore also the 
corresponding Ermakov invariant) are identical to those of our logarithmic NLSE (69). 
There is a small point that cannot be explained sufficiently in Hasse's approach, namely, 
a contribution to the energy originating from the non- vanishing mean value of Wn a s, 

(w Has ) / o. 

Comparison of Whcis^wp(x, t) with our logarithmic term Wsch^wp(x, t) shows that 

WsCH = Wuas - (W Has ), (79) 

so, in our case, the non-vanishing mean value does not occur. The relationship between 
this approach using the "friction potential" and our logarithmic ansatz is connected with 
a unitary transformation and is discussed in more detail in [34]. 

In the following, the logarithmic approach will be considered explicitly since the rele- 
vant part for the CS is essentially the complex Riccati equation (70) that is identical to 
Hasse's approach. But it is more straightforward to bridge the gap between the afore- 
mentioned explicitly TD Hamiltonians based on non-canonical modifications of classical 
mechanics and the nonlinear modifications of the SE by using the complex logarithmic 
term. 



3.5 Connection between the canonical and nonlinear approaches 

To establish the connection between the explicitly TD Caldirola-Kanai approach and the 
logarithmic NLSE (69) we refer to Schrodinger's first communication on wave mechanics 
[1] where he starts from the Hamilton-Jacobi equation 

+ =0 (80) 

with the action function S and the momentum p = -JLS. He introduced the wave function 
ty(x,t) via S c = jln\l/, where the subscript c (added by us) indicates that this action is a 
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complex quantity, since * is, in general, a complex function (a fact that Schrodinger did 
not like at all in the beginning |37J). Via a variational ansatz, Schrodinger arrived at the 
Hamiltonian operator Hl = — + V (x). 

We now reverse Schrodinger's procedure, starting with Eq.(69) (divided by *, which 
causes no problems for Gaussian WPs) and, using the definition of S c , arrive at 

^ + 1^)S C + H = - 1 (S C }. (81) 

This is, of course, as little rigorous as Schrodinger's first attempt was. However, it follows 
his idea of connecting the classical Hamilton-Jacobi theory with a wave (mechanical) 
equation. The purely TD term —j(S c ) is necessary mainly for normalization purposes 
(can therefore be absorbed by the normalization coefficient) and shall be neglected in the 
following. 

Multiplying the remaining Eq.(81) by e 7 * and using the definitions 

S c = e^S c and H = e^H (82) 
it can be rewritten as Hamilton-Jacobi equation 

^S c + H = 0. (83) 

From the definition of the action function, it follows that the wave function ty{x,t) in 
the transformed (canonical) system is connected with the wave function ty(x,t) in the 
physical system via the non-unitary relation 

In* = e 7 *ln*. (84) 

Consequently, the (complex) momenta in the two systems are connected via 

p c = -—In* = e 7 '- — ln*^ = e 7 V c , (85) 

l OX l ox 

which is equivalent to the connection between the canonical and the kinetic momentum 
in the Caldirola-Kanai approach. The non-canonical connection between the classical 
variables (x,p) and (x,p) corresponds to the non-unitary transformation between * and 
*. 

Note: Although * and * depend explicitly on the same variables, x and t, the two 
wave functions are analytically different functions of x and t and have different physical 
meanings due to the non-unitary transformation (84). Ignoring this fact leads to the 
apparent unphysical results like violation of the uncertainty principle (for further details, 
see P^). 

Expressing H in terms of the canonical momentum p c and following Schrodinger's 
quantization procedure, finally yields the modified SE (53) of the Caldirola-Kanai ap- 
proach. 
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Because the WP solution of Eq.(53), ^/ck(x, t), is related to the WP solution ^nl(x, t) 
of the NLSE (69) via (84), this means particularly for the corresponding Riccati and 
Ermakov equations that 



(^ 



= e 



it 



CK 



2ft 
—y 

m 



NL 



and 



olck 



(86) 



(87) 



Taking this into account, Eqs.(55) and (56) turn into Eqs.(70) and (72) and the Ermakov 
invariant Iqk turns exactly into I^l- 



3.6 Uncertainties, bifurcations, creation and annihilation oper- 
ators and coherent states for dissipative systems 

With the help of the WP solution of Eq.(69) the uncertainties of position and momentum, 
their correlation and product as well as their contribution to the energy of the system can 
be expressed in terms of a^L and a^L or (^Ui) NL and (^Ur) NL , respectively, where the 
latter ones are more general since this form can be used identically for the corresponding 
quantities in the case without dissipation (therefore the subscript NL is omitted). One 
obtains 

ho ft ( 2ft 



2m 



-a 



NL 



2m 



m 



-yi 



(88) 



(P 2 



'NL 



{p-NL - |«Vl) 2 + 7T- 



¥(^) _1 kf^) 2 +&) 2 l, (89) 



([x,p\+) NL = ft (a NL a NL - t^tvl) = ft (~y j 



2ft V 1 (2ft 



m 



-Vr , 



(90) 



U NL = i^) NL (p% L = i\l+ {a NL a NL - la% L ) 2 ] = J [l + (f yj)' 1 {^y R f] (91) 



E 



NL 



2m 



/-~2\ 2 i-~- f l\ ft 

(p )nl + (a; }nl = j 



OiNL~-OiNLl H 5 rUJ a NL 



a 



NL 



ft (2ft 



= 7 ~ Vl 
4 \ m 



-i 



2ft Y (2ft 

—Vr + —Vi J + W 
ml \m 



(92) 



At this point we wish to recall the solution of the Riccati equation (4) via transformation 
into a Bernoulli equation (9) with the help of a particular solution of the Riccati equation, 
where this solution (^y) occurred in the Bernoulli equation as parameter A = 2 (^y)- 
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Including the dissipative term in Riccati equation (70), this parameter (already in the 
simples cases for constant (—y)) now takes the form 



.4 



1/2 



(93) 



So, already for the free motion [oj = 0), the two values A± = ±7 are possible, leading 
to two different solutions with different (TD) uncertainties and contributions to the en- 
ergy, Enl±, as can be seen from Eqs. (88-92). This kind of Hopf bifurcation also occurs 
for oj 7^ and can be expected due to the quadratic nonlinearity of the Riccati and 
Bernoulli equations; however, it is not immediately obvious when considering the equiv- 
alent Ermakov equations. Bifurcations of this kind are common in classical nonlinear 
dynamics [38], but usually in connection with real evolution equations. In our quantum 
mechanical context, the bifurcation does not show up in the conservative case because 
either the two solutions are degenerate (V = : A = ±0) or one of the two leads to 
unphysical results (HO : A = —i2uo). Taking into account the dissipative environment 
changes this situation qualitatively. Different solutions of the Riccati equation (70) also 
entail different creation/annihilation operators for the respective system. 

In the conservative case we found that the creation and annihilation operators are 
related to the Ermakov invariant via 



h 

m 



a + (t)a(t) + 



1 



h 



or 



-zz 



m 



(94) 



i.e., they could also be obtained by factorization of /imp- 
leads to 



Trying to do the same with Inl 



znl 




— ^— J + i (r]a NL - (a NL - ^-a NL ) 77 
a NL l \ \ 2 



7/ 



I { '/'> vi - ( a N L - TjOiNL j >7 



(95) 



or 



a NL (t) 



NL 



(t) 




Pop ( 2h 

m \m 



x 



NL 



(96) 



Pop 

m 



2h t 
—y 

m 



x 



NL 



From (96a,b) it becomes obvious that, apart from the exponential factor e 7 */ 2 , the form 
of the creation and annihilation operators is identical to the ones defined in (21a,b). Only 
the solutions of Riccati equation (4) must be replaced by the ones of Eq.(70) belonging 
to the logarithmic NLSE (69). 

Still, the meaning of the exponential factor e 7 */ 2 needs to be explained. Since Inl 
is a constant of motion for the dissipative system, like H exp on the canonical level, but 
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the nonlinear Hamiltonian and its mean value are not, the invariant should correspond to 
the canonical level that is related to the physical one via the non-unitary transformation 
ln\l/ = e 7 *ln\l/. Therefore, the factor e 7 * in I ml should be omitted if it is used for the 
definition of the creation and annihilation operators that supply the wave functions and 
CSs which are solutions of the Hamiltonian operator on the physical level. 

So, finally the generalized creation and annihilation operators for the dissipative sys- 
tems can also be written in the form 



a(t) 
a + (t) 



m 



2h \mdx \m 





in i h d 

2h \ mdx 



2h 



m 



x 



(97) 



The corresponding CSs can again be obtained via application of a(t) as eigenstates of this 
operator or via the displacement operator, leading to 



--\z\ 2 



n=0 



z n (a H 



nl 



-|o>, 



(98) 



with 




2h 



m 



Vi + 



2h 



m 



-Vr V 



where the CSs \z) are identical to the exact WP solutions ^wp of the NLSEs discussed 
in Sections 3.3 and 3.4 and related to the canonical WPs of Sections 3.1 and 3.2 via the 
non-unitary transformation given in Section 3.5 (for further details, see also [43]). 



4 Conclusions 

A generalization of CSs and the creation/annihilation operators used to construct them 
is necessary if they represent WPs with time-dependent width corresponding to TD po- 
sition uncertainties. This situation not only occurs for HOs with TD frequency u = uj(t), 
but already for oj = u = constant if the initial state is not the ground state of the 
harmonic oscillator. In these cases, the CS is not a minimum uncertainty CS (since po- 
sition/momentum correlations are involved) but still fulfills the Schrodinger-Robertson 
uncertainty relation. Using a Gaussian WP ansatz, from the corresponding TDSE the 
equation of motion for the WP width can be obtained via a nonlinear Riccati equation 
for a complex variable (—y(t)). This equation can be transformed into an Ermakov equa- 
tion for a quantity a(t) which is directly proportional to the WP width. This equation, 
together with the Newtonian equation that describes the classical motion of the WP max- 
imum, allows for the definition of the dynamical invariant 7^. Instead of considering the 
Ermakov equation, we concentrate on the complex Riccati equation and its formal analyt- 
ical solutions which not only depend on the physical system under consideration, but also 
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qualitatively on the (complex) initial conditions. These solutions have close formal sim- 
ilarities with general superpotentials leading to isospectral potentials in supersymmetric 
quantum mechanics. 

The generalized creation/annihilation operators can be expressed in terms of the com- 
plex Riccati variable using the relation between its imaginary part and uo occurring in the 
standard form, i.e., (—yA = 0Jq, and replacing w in a and a + by the complex quantity ac- 
cordingly. The resulting operators are in agreement with those obtained via factorization 
of the operator corresponding to II. 

Using the definition of the CS as eigenstate of the annihilation operator, it is possible 
to determine the complex eigenvalue z(t) in terms of the classical quantities i] and rj and 
of real and imaginary parts of (— y). The CS in position space obtained in this way is 
identical to the Gaussian WP that is a solution of the corresponding TDSE. 

In addition, the definition of the CS via the displacement operator, expressed using 
the generalized a(t) and a + (t), provides the same Gaussian WP. Finally, the mean value of 
position and momentum, as well as their uncertainties, are calculated using the generalized 
CSs and it is proven that the Schrodinger-Robertson uncertainty condition is fulfilled. 

In the second part, the formalism is extended to include dissipative environmental 
effects. For this purpose descriptions of the dissipative systems in terms of modified SEs 
are considered. Specifically, two approaches are discussed based on classical Hamiltonians 
obtained via non-canonical transformations leading to modified, but linear, TDSEs on a 
formal canonical level and are compared with two approaches using nonlinear modification 
of the TDSE. In all cases, Gaussian WP solutions again exist and the relevant complex 
Riccati equation is the same for both nonlinear approaches. The Riccati equations of the 
canonical approaches are related via a unitary transformation of the corresponding WPs 
and can be transformed into the one of the nonlinear approaches taking into account the 
non-unitary transformation (84) between the canonical and the physical level (where the 
NLSEs are valid). 

In the dissipative case the mean value of position, and hence the WP maximum, obeys 
a Newtonian equation with additional friction term linear proportional to velocity. The 
position and momentum uncertainties, and their contributions to the energy of the CS, 
are also determined. When expressed in terms of the auxiliary Ermakov variable Q-nl (t) 
they have a different form from the case without dissipation, particularly because the 
relation between (^Vr) and a, a is different. However, expressed in terms of real and 
imaginary parts of the form is identical to the one in the conservative case. The 

same applies to the creation/annihilation operators that can be obtained in the same way 
as in part one. They also agree with those obtained via factorization of the Ermakov 
invariant corresponding to the dissipative systems (if the difference between canonical 
and physical levels is taken into account). 

Concluding, one can say that creation and annihilation operators, corresponding CSs 
and quantum uncertainties for time-dependent conservative and dissipative Hamiltonians 
can all be expressed in the same form if the complex variable (— y\ of the Riccati equa- 
tion is applied. Due to the nonlinearity, for the same physical system different CSs with 
different physical properties can exist according to the respective solution of the Riccati 
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equation. Furthermore, these solutions can also be very sensitive to the (complex) ini- 
tial conditions. In this context, formal similarities and comparison with supersymmetric 
quantum mechanics and the CSs of the nonlinear algebras involved [391 - H2] might be useful 
and will be investigated. 
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